Resonant states and pseudo-spin symmetry in the Dirac Morse potential 
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The complex scaling method is applied to study the resonances of a Dirac particle in a Morse 
potential. The applicability of the method is demonstrated with the results compared with the 
available data. It is shown that the present calculations in the nonrelativistic limit are in excellent 
agreement with the nonrelativistic calculations. Further, the dependence of the resonant parameters 
on the shape of the potential is checked with the sensitivity to the potential parameters analyzed. 
By comparing the energies and widths of the pseudospin doublets, well pseudospin symmetry is 
discovered in the present model. The relationship between the pseudospin symmetry and the shape 
of the potential is investigated by changing the Morse potential shaped by the dissociation energy, 
the equilibrium intermolecular distance, and the positive number controlling the decay length of the 
potential. 
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I. INTRODUCTION 



O !■ It is well known that the resonances are the most striking phenomenon in the whole range of scattering experiments, 
and appear widely in atomic, molecular, nuclear physics and in chemical reactions. One of the best-known examples 
is the spontaneous positron production in uranium nuclei collisions, where some harp resonance peaks have been 
observed [l|. 

p~ ! To explore the resonances, ones have developed many techniques, which include the R-matrix theory Q, K- matrix 
i . theory 0], the scattering phase shift method [J|, and several bound-state-like methods, such as the real stabilization 
method (RSM) Q, the analytic continuation in the coupling constant (ACCC) method Q and the complex scaling 
CN , method (CSM) [3, Q. Due to the simplicity in the calculations, these bound-state-like methods have been widely 
employed to study resonances in quantum many-body systems. For examples, much effort has been made in order to 
calculate more efficiently resonance parameters with the RSM [9j-[l2j. Combined with the cluster model, the ACCC 
■ approach has been used to calculate the resonances in some light nuclei (13 , Il4j . In combination with the relativistic 
\ mean field (RMF) theory [15 -Tl71| , t he RMF-ACCC has presented a good description for the structure of resonant 
levels in some realistic nuclei 18l Il9j. Moreover, the CSM have been successfully used to study resonances in atomic 
and molecular systems [2(J [H | and atomic nuclei {22M2t| . 
. Recently, the Morse potential [13] attracts additional attentions for its success in describing the motion of atomic 
y—( ' and molecular as well as nuclear systems. The vibration of a diatomic molecule was well depicted by the Morse 
potential |28| . The Morse potential is also very useful in expressing the properties of nuclei [29l - [31| . Other applications 
of Morse potential include the description of the interactions between two atoms in a diatomic molecule [32H35j , the 
interatomic potential of crystalline solids [36| . and the potential for the adsorption of a molecule or atom by a solid 
surface [13] ■ This list is far from comprehensive and new applications continue to be found as well. 

Due to the wide applicability, it is interesting to probe the resonances in a Morse potential. Several methods 
have been developed to study the resonant states of a nonrelativistic particle in a Morse potential. Nasser et al. [H[ 
have adopted the J-matrix method to examine not only the bound states but also resonances associated with the 
rotating Morse potential model. Satpathy et al. [391442T ] have employed the long-ranged Morse potential between two 
carbon nuclei in order to obtain many vibrational states which can be compared with observed resonance levels. Kato 
and Abe [43| have applied the Morse potential model to the 12 C- 12 C system as a typical example of the molecular 
resonances observed in many lighter heavy-ion collisions. Jarukanont et al. 4J] have applied the Floquet theory in 
combination with the exterior complex scaling to obtain the energies and the distributions of probability for the 
quasibound states of a driven Morse potential. 

In addition to these nonrelativistic studies, much attention has been paid to the problem of the relativistic Morse 



potential. In 2001, the s-wave Dirac- Morse problem was formulated and solved exactly [45[. As there are no analytical 



solutions to the Dirac equation with k^O for the rotational Morse potential, some approximations are employed to 
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present the numerical or quasi-analytical solutions [46|,|47(. To recognize the rclativistic effects of the Morse potential, 
Dirac equation has been solved for attractive scalar and repulsive vector Morse potentials under pseudospin symmetry 
by the Pekeris approximation [48], |49| . Other interesting work includes the solution of the relativistic Morse potential 
problem under the condition of pseudospin symmetry |50( ■ 

These relativistic studies on Morse potential associated with an interesting concept, pseudospin symmetry (PSS). 
The concept was first introduced in the field of nuclear physics more than 40 years ago [5lL |52j . Its meaning is that the 
single particle states with quantum numbers (n, I, j = 1 + 1/2) and (n — 1, 1 + 2, j = 1 + 3/2) are nearly degenerate. Until 
1977, this symmetry was known as a relativistic symmetry and the condition of exact PSS was proposed in Ref. [53j . 
Afterwards, the condition of exact PSS was extended to a more general case (54|, and found to be approximately 
satisfied in exotic nuclei with highly diffuse potentials [55| . The dynamical character of the symmetry was recognized 
in Ref. [56j and the spin symmetry in the antinucleon spectrum was indicated in Ref. (57j . Several recent progresses 
include the supersymmetric description for deformed nuclei [58[, the nonperturbation nature of PSS [53 . |6G |. the 



relativistic effect of PSS [61j, the influence of tensor interaction on PSS [62, [63|, the role of central Coulomb potential 
at PSS [iij], and the spin and pseudospin symmetries of the Dirac equation with confining central potentials [65| . In 
Ref. [f36j], the similarity renormalization group is used as the critical tool for understanding the origin of the PSS and 



the cause why the PSS becomes better for the levels closer to the continuum is disclosed in a quantitative way [67 1. 
In combination with the supersymmetric quantum mechanics, the origin of the PSS and its breaking mechanism are 
detailedly analyzed at the nonrelativistic limit 1681 . This symmetry is also checked in the resonant states with similar 
features to bound states indicated in Ref. (69l . l70l |. where it was indicated that in conditions of pseudospin symmetry 
the same pseudospin quantum numbers will be conserved and the pseudospin doublets would have the same energy 
and width. As the interesting phenomenon also appears in the resonant states, in this work we will apply the CSM 
to study the resonance states of a Dirac particle populating in the Morse potential, and survey the existence of PSS 
for this system. 



II. THEORY 



The Dirac equation of describing a particle moving in the Morse potential V(r) is 
with 

V(r) = V a e~ {r ~ ro)a (2 - e -( r - r °M , (2) 

where k = +(j + 1/2) for j = 1 + 1/2, Vo is the strength of the potential (or the dissociation energy in the context of 
diatomic molecules), rg is the equilibrium intermolecular distance, and a is the positive number controlling the decay 
length of the potential. For Vb < 0, this potential has a minimum of Vo at r = and it is called the regular Morse 
potential. For Vo > , it has a maximum of Vb there and it is called the inverted Morse potential, and then r$ is the 
location of the top of the barrier. 

The starting point of the CSM is a transformation of the Hamiltonian H , which is operated by defining an unbounded 
nonunitary scaling operator U(9) with a real parameter 9. Similar to Ref. [25L [7lj . the complex scaling operator takes 
the form 

/ e l9§ \ 

where S = r-§-. Then, the Dirac Hamiltonian in Eq. ([TJ becomes a complex scaling form: 



Ho-U{V)HU [0) - I e _ l6 ( d + l±K\ y re i6\ _ M ■ (4) 



The corresponding complex scaled equation is 

Hgtp e {r) = seifje (r) , (5) 

where ipg = U (0) ip (r) is the complex scaled wave function. According to the Aguilar-Balslcv-Combes (ABC) theorem 
[72j, the energy for the bound states is independent of 9, and the real part E r and the imaginary Ei part of the energy 
Eg for the resonant states are independently of 9 as well. T = —2Ei represents the width for the resonant states. 
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In this work, Eq.([5|) is solved by expanding the large and small components of the Dirac spinors / (r) and g(r) in 
terms of the radial functions R n i (r) of a spherical harmonic oscillator (HO) potential, i.e., 

fe (r) =J2fn (0)R n i (r) ; ge (r) = £ gn {6)R h l (r) , (6) 

n— 1 n— 1 

where the orbital angular momenta / and I have the same meaning as that in Ref. [73| . The upper limits n max and 
ii max in Eq.(6) are radial quantum numbers determined by the corresponding major shell quantum numbers -ZV max 
= 2 (n ma x — 1) + 1 and -/V max = 2 (n max — 1) + Z. In the present calculations, the maximum of the major shell quantum 
number N = 200 is chosen to ensure the calculated result with the accuracy better than the order of 10 -4 fm -2 or 
10~ 4 a.u. for the energies and widths. The details can be referred to the literatures (7ll l73j. 

Inserting Eq. ([6]) into Eq. ([5]) and using the orthogonality of wave functions R n i , one arrives at a symmetric matrix 
diagonalization problem of the dimension n, max + n max , i.e., 



V ni n> + M ■ 5 n>n i B n fii \ ( f n r \ _ / /„ 

fls.»' Vfl «/ - M • 5* ft/ / I g w I ~ ° e I ,9ft 



(7) 



The matrix elements V n ,n' and ^n,^ 7 are given by 

V n , n , = / r 2 dr R nl {r) [V{re w )} R n u{r), (8) 



r dr 



d 1 + k 



dr r 



Rn>i(r) 



(9) 



The integral in Eq.(8) is calculated with the Gauss quadrature approximation [20l. l71|. The matrix elements B n ^ 
can be further simplified as 



flfi,n/ — 



{^/fl + l + Tj28n,n' + V^^.n'-l) , K < 



(10) 



k^i<r f -y/n + Z - l/2^ft,„' + \/n - 15ft,„/+iJ , k > 
With the matrix elements Vnn' an( i A« fi'; the solutions of the Dirac equation are obtained by diagonalizing the matrix 



He 



(ii) 



The eigenvalues of Hg representing bound or resonant states do not change with 8, while eigenvalues representing 
continuous spectrum rotate. 



III. RESULT AND DISCUSSION 



With the theoretical formalism represented in the previous section, we explore the resonant states for a particle 
populating in a Morse potential. For comparison with the results from the J-matrix approach in Ref. [38| . the same 
values for the parameters of the Morse potential are employed. The illustrated results are plotted in Fig.l with the 
complex rotation angle = 70° for the states with k = — 1. It should be mentioned that the energy displayed in Fig.l 
is multiplied by a factor 2M/h 2 in order to compare with Ref. [IU. The same as Fig.l, the energy and the depth of 
potential in the following Table I and Figs. 2 and 3 are multiplied by a factor 2M/H 2 . All the eigenvalues of Hg, which 
correspond to the bound, resonant, and continuum states, are respectively labeled as open boxes, open circles, and 
solid circles. From Fig.l, one sees clearly that the eigenvalues of Hg fall into three regions: the bound states populate 
on the negative energy axis, while the continuous spectrum of Hg rotates clockwise with the angle 29, and resonances 
in the lower half of the complex energy plane located in the sector bound by the new rotated cut line and the positive 
energy axis get exposed and become isolated. In the present calculations, as the finite basis is used, the continuous 
spectrum of Hg consists of a string of points. 

The energies for the resonant states in Fig.l are listed in Table |U The rclativistic results are shown in the first 
column. The second column displays the data in the nonrclativistic limit, which are obtained with an approximation 
by increasing the value of the speed of light up to 100 times larger than its actual value in the relativistic CSM 
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calculations [20(. In order to compare with the other methods, some data from the nonrclativistic calculations are 
exhibited in the third column by J-matrix method [38j j and the fourth column by S- matrix method 74 1- From Table HI 
it can be seen that the deviation between the relativistic result and that in the nonrclativistic limit is observable, 
which implies the relativistic effect can not be completely ignored in the present model. The present calculations 
in the nonrelativistic limit agree with the nonrelativistic calculations quite well. The deviation between our results 
and the data by J-matrix or S-matrix method is very small. Whether for a relativistic or nonrelativistic system, the 
present method has provided a good description for the resonant states. Hence, we can look into all the resonant 
states for a Dirac particle in a Morse potential by this method. 

To survey the resonances in a wider range, it is necessary to analyze the dependence of the resonant parameters on 
the shape of the potential, which can help us to recognize the resonances for more realistic systems. The variations 
of energy and width with the parameters of the potential are shown in Fig. 2, where the data associated with the 
dissociation energy, the equilibrium intermolecular distance, and the positive number controlling the decay length of 
the potential are respectively shown in the panels (a), (b), (c) of Fig. 2. Starting from Vb = 1.5 fm~ 2 , the energy of 
the resonant states increases with the increasing of Vq. Further increasing Vb, the energy level changes the direction 
of evolution for the lsi/2, followed by the 2s±/2, and so on. Simultaneously, some new resonant levels are exposed. 
For the width, its evolution with Vb is consistent with that for the energy. Starting from Vb = 1.5 fm~ 2 , the width 
reduces with the increasing Vb for the several lower levels. For the several upper levels, the width increases with the 
increasing of Vo- As n increases, the width is sequentially changed in the direction of the evolution with Vq- The 
phenomenon can be explained by the potential energy curves plotted in Fig. 3. With the increasing of Vb, the barrier 
becoming higher and the well becoming deeper arc seen in Fig. 3(a) for the Morse potential. It is the cause that the 
resonant parameters vary with Vb in accordance with the foregoing. With the increasing of n, the evolution of the 
resonant energy with ro changes from a decreasing to an increasing with the increasing tq. The level located on the 
middle, the energy of the state (4s!/ 2 ) only has a slight change. The energy for the levels lower than 4si/2 decreases 
with the increasing ro, while the energy increases with the increasing ro for the energy level populating in higher than 
4si/2j which is associated with the variation of the potential. From the potential energy curves in Fig. 3(b), one can 
see that both the barrier and the well become wider with the increasing of ro, which leads a weakening of penetration 
for a particle in the Morse potential. Thus, the width becomes systemically smaller with the increasing ro. Fig. 2(c) 
shows the resonant parameters varying with a, a similar change is seen for the energy with ro in Fig. 2(b) and for 
the width with Vb in Fig. 2(a). These can be explained from the changes in the potential field. The potential energy 
curves in Fig. 3(c) display that the potential well becomes deeper and the range of potential barrier becomes narrower 
with the increase of a, which is the cause of the resonant parameters with a. 

To further examine the present CSM calculations, the resonant states with k = 2 is investigated with the data 
exhibited in Tabic [TT1 For comparisons with the data in Ref. |38| . the atomic units h = m = 1 is adopted here. Similar 
to Table HI the result in the nonrelativistic limit is obtained according to the method in Ref. [2(| . From Table UH it 
can be seen our results are in excellent agreement with the data in Ref. [H[ , which are obtained by J-matrix method 
for a nonrelativistic system. These show the present calculations is fully correct, and we can explore the resonant 
states of a relativistic particle by this method. 

With the resonance states obtained by the CSM, the position of pseudospin doublets in the complex energy surface 
is shown in Fig. 4, where the pseudospin doublets k = — 1 and 2, K = —2 and 3, K = —3 and 4, as well as K = —4 and 5 
are respectively displayed in Fig. 4(a), (b), (c), and (d). It is clear that PSS is well preserved for all the doublets. For 
the pseudospin doublet with the same orbital angular momentum, the quality of PSS improves with the increasing of 
the radial quantum number. 

To disclose the relationship of the PSS in the resonant states and the shape of the potential, we analyze the variation 
of the pseudospin splittings with the parameters in the Morse potential. Keeping ro and a fixed, we vary Vb to see how 
the energy and width of the pseudospin doublets are sensitive to the dissociation energy. This dependence is shown 
in Fig. 5(a) for the energy and width splitting. It is clear that the energy splitting of the pseudospin doublets increase 
with the increasing of Vo, while a inverse trend is found for the variation of width. This can be easily understood as 
higher barrier and deeper potential well are accompanied with the increase of Vq. Similarly, the energy and width 
splittings varying with ro is displayed in Fig. 5(b) with the other parameters fixed. From the left (ro = 0.5), as the 
radius ro increases, the energy and width splittings decrease, which can be understood by the dependence of dV(r)/dr 
on r. The variation of the pseudospin splittings with a is presented in Fig. 5(c). For all the pseudospin doublets, 
the tendency for the change of the pseudospin energy and width splittings with a is the same. With a increases, 
the splitting increases. This could be expected because the potential well becomes deeper and widened with the 
increasing of a, hence the increase of the derivative of V(r). In the range of the potential parameters considered here, 
the deviations between pseudospin doublets is within 1.17 a.u. for E r and 1.04 a.u. for V, which indicates the PSS is 
well preserved. 
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IV. SUMMARY 

In summary, the resonant states of Morse potential are investigated by using the complex scaling method in the 
rclativistic framework. The calculated results are compared with the available data in references and satisfactory 
agreements are found. Further, the dependence of the resonant parameters on the shape of the potential is checked 
and a unified result is found for the evolution of the resonant parameters with the potential parameters, which can 
be well explained from the change of potential energy curves. With these resonant states obtained, the PSS in the 
resonant states is investigated and well PSS is found in the present model. The pseudospin splitting is shown in 
correlation with the Morse potential shaped by the dissociation energy V , the equilibrium intermolecular distance r , 
and the positive number controlling the decay length of the potential a, in which three Morse potential parameters 
Vq, a and ro are found to play the important roles in the splittings of energy and width of pseudospin doublets. 
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FIG. 1: (Color online) Position in complex energy surface for the states with « = —1, where the bound, the resonant, and the 
continuum are respectively labeled as open boxes, open circles, and solid circles. The parameters in the Morse potential are 
adopted as Vo = 6 fm~ 2 , ro = 4 fm, and a = 0.3 fm -1 . The mass of particle M = 0.5 fm -1 and the rotation angle 9 = 70° are 
adopted in the present calculations. 




2 4 6 8 10 2.8 3.4 4.0 4.6 5.2 0.18 0.24 0.30 0.36 0.42 



V [fm" 2 ] r o [fm] a [fm 1 ] 

FIG. 2: (Color online) Energies and widths as a function of every potential parameter for the resonant states shown in Fig. 1. 
The data corresponding to the variables Vo, ro, and a are respectively displayed in the subfigures (a), (b), and (c) with the 
other parameters fixed to ro = 4 fm, a = 0.3 fm -1 in (a), Vo = 6 fm -2 , a = 0.3 fm -1 in (b), and Vo = 6 fm -2 , ro = 4 fm in 
(c). The mass of particle is adopted as M = 0.5 fm -1 . 
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FIG. 3: (Color online) The Morse potential model for (a) various strength of the potential Vo = 1,6, 10 fm -2 along with 
ro — 4 fm and a = 0.3 fm -1 , (b) various equilibrium intermolecular distance ro = 3,4,5 fm along with Vo — 6 fm~ 2 and 
a — 0.3 fm _1 , and (c) various positive number a = 0.2, 0.3, 0.4 fm -1 along with Vo = 6 fm -2 and r = 4.0 fm. 



I-' 

-5 
-10 

-15 I- 



E ..20 
-5 



-10 - 



-15 - 



-20 - 



1 1 1 I 1 1 

• K=-l 

O K=2 



•O 2s lo ,1d„ 



O 3s , 2d 

^ 1/2' 3/2 



O 4s , 3d 

1/2' 3/2 

O 5s , 4d 

^» 1/2' 3/2 



r> 6s , 5d 

^ 1/2' 3/2 

© 7s , 4d 

1/2' 3/2 



(a) :: 



O 6 P3,2' 5f 5 , 



• K=-3 

O K=4 



•o 2d , 1q 

"SQ' »7/2 

^> 3d , 2q 

5/2' B 7/2 



^> 4d , 3q 

5/2' a 7/2 

^ 5d 5,2' 4 9 7 , 2 
^ 6d 5,2' 5 9 7 ,2 (C) 



_L 



_L 



_L 



1 " 1 I 1 

• K=-2 

O K=3 



T 



I 1 1 1 LI 



• 1 p 3 , 



•O 2p 3(2 , 1f 5(2 



*> 3p 3(2> 2f 



«> 4 P3,2'3f 5 , 



O 5 P 3 /2' 4f 5,2 



3/2' 512 
P 7 P3/2'. 6f 5,2 



(b) : 



j i i i i i i i i i i i i i i i i i i ~| i i i i i i i i i i i i i i i i i i r 



• K=-4 

X~ o k=5 



• 1f 7,2 

•b 2f 7 ,2' 1h Si 

3f 7 ,„2h a „ 



*> 4f 7 ,2' 3h 9 /2 



^ 5f 7 ,2' 4 \ 



(d) : 



10 



15 



10 



15 



20 



FIG. 4: (Color online) Pseudospin doublets in the complex energy surface for the Morse potential with the parameters Vo = 10, 
ro = 1, and a — 0.5 in the atomic units h = m = 1. 
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TABLE I: Energies in the present calculations for the states with « = — 1 in comparison with the data in Refs. |38|, UM- Here 
the parameters of Morse potential are adopted as Vb = 6 fm -2 , ro = 4 fm, and a = 0.3 fm -1 . The mass of particle is adopted as 
M = 0.5 fm -1 . For comparison with Refs. [H, [74|]), the depth of potential displayed and the calculated energies are multiplied 
by a factor 2M/H 2 . 



E r + iEj( fm~ 2 ) E r + iEj( fm~ 2 ) E r + iEj( fm~ 2 ) E r + iEj( fm~ J ) 



Relativistic 



Nonrelativistic limit 



Nasser 



Rawitscher [741 



-8.1096 



-8.1089 



1.1745 
5.6229 
6.8906 
7.3194 
7.1143 
6.3679 
5.1514 - 
3.5200 - 
1.5173 - 
-0.8212 
-3.4656 
-6.3907 
-9.5750 
-12.9998 
-16.6491 
-20.5088 



-i 0.0002 
-i 0.0349 
-i 1.3170 
-i 3.5858 
-i 6.0693 
-i 8.5999 
i 11.0980 
i 13.5209 
i 15.8439 
-i 18.0518 
-i 20.1350 
-i 22.0869 
-i 23.9031 
-i 25.5801 
-i 27.1155 
-i 28.5077 



1.1779 
5.6212 
6.9041 
7.3375 
7.0545 
6.3820 
5.1806 - 
3.4649 - 
1.5262 - 
-0.8167 
-3.4773 
-6.3930 
-9.5671 
-12.9859 
-16.6258 
-20.4731 



■i 0.0011 
■i 0.0303 
■i 1.3103 
■i 3.6189 
■i 6.0853 
■i 8.5374 
i 11.1429 
i 13.5147 
i 15.8071 
1 18.0490 
i 20.1177 
i 22.0537 
i 23.8658 
-i 25.5370 
-i 27.0640 
-i 28.4487 



1.1778 
5.6252 
6.8911 
7.3182 
7.1111 
6.3627 
5.1446 - 
3.5123 - 
1.5095 - 
-0.8278 
-3.4697 
-6.3907 
-9.5691 
-12.9861 
-16.6255 



.1090 

-i 2.01E-13 
i 0.0351 
i 1.3194 
i 3.5887 
i 6.0715 
i 8.6005 
11.0960 
13.5151 
15.8334 



-8.1090 
1.1783 
5.6252 -i 0.0351 



i 
i 
i 

-i 18.0358 
-i 20.1128 
-i 22.0579 
-i 23.8667 
-i 25.5363 
-i 27.0642 



TABLE II: The same as TABLE I, but for the states with k = 2 in comparison with the data in Ref. [3^|. Here the parameters 
of Morse potential are adopted as Vb = 10, ro = 1, and a — 2.0 in the atomic units h — m — 1. 



E r + iE t 


E r + iE t 


E r + iEi 


Relativistic 


Nonrelativistic limit 


Nasser [38] 


-30.7047 


-30.4136 


-30.4139 


10.8020 -i 0.2822 


10.9262 -i 0.3026 


10.9260 -i 0.3027 


17.1419 -i 12.3689 


17.1244 -i 12.5031 


17.1240 -i 12.5027 


11.1795 -i 32.0868 


11.0511 -i 32.1914 


11.0521 -i 32.1906 


-4.8377 -i 52.5443 


-5.0383 -i 52.5395 


-5.0376-i 52.5407 


-29.3283 -i 72.2381 


-29.5208 -i 72.0565 
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FIG. 5: (Color online) Pseudospin energy and width splittings as a function of every potential parameter for the resonant states 
shown in Fig. 4(a). The data corresponding to the variables Vo, ro, and a are respectively displayed in the sub-figures (a), (b), 
and (c) with the other parameters fixed to ro = 1, a = 0.5 in (a), Vo = 10, a = 0.5 in (b), and Vo = 10, ro = 1 in (c) in the 
atomic units H = m = 1. 



